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ABSTRACT 

We have simulated the 3D emergence and interaction of two twisted flux tubes, 
which rise from the interior into the outer atmosphere of the Sun. We present 
evidence for the multiple formation and eruption of flux ropes inside the emerging 
flux systems and hot arcade-like structures in between them. Their formation 
is due to internal reconnection, occurring between oppositely directed, highly 
stretched and sheared fieldlines at photospheric heights. Most of the eruptions 
escape into the corona, but some are confined and fade away without leaving the 
low atmosphere. As these flux ropes erupt, new reconnected fieldlines accumulate 
around the main axis of the initial magnetic flux systems. We also show the 
complex 3D fieldline geometry and the structure of the multiple current sheets, 
which form as a result of the reconnection between the emerging flux systems. 
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1. Introduction 



Partial eruption of a flux rope can be produced by favourable, imposed shearing motions 
at the photosphere. This shearing can occur as sunspots with opposite polarity move apart 
along the neutral line of an active region. Then, as the rising magnetic system expands, the 
sheared magnetic fieldlines are stretched and undergo internal reconnection, at low heights. 



which leads to partial eruption of magnetic flux (Manchester et al. 2004). This is consistent 



with the 'tether-cutting' mechanism (Sturrock 1989; Moore et al. 2001), in which, effective 



reconnection of magnetic fieldlines at the core of a magnetic bipole releases the energy stored 
in the sheared field, resulting in an eruption. An important factor in determining whether 
an emerging flux rope will erupt is the amount of twist of the fieldlines around the central 



axis of the tube (Murray et al. 2006). If the winding of the fieldlines is high enough, a kink 



instability occurs that eventually leads to either a confined or an ejective eruption (Torok & 



Kliem 2005; Gibson & Fan 2006) 
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Observations ( Sterling &: Moor"e||2004 2005 Jiang et al.||2007 ) have shown active-region 
filament eruptions when new flux emerges near the location of a preexisting emerging flux 
region. It is likely that the new emergence is important in the eruption onset. Several nu- 
merical models, requiring multipolar magnetic conflgurations for eruptions, have reproduced 



the above and similar-like observations. In the breakout model (Antiochos 1998), a sheared 
bipole rises into an overlying bipole and 'external' reconnection occurs at a magnetic null 
point at their interface. When the reconnection becomes fast enough the rising bipole breaks 
through the ambient fleld and erupts. 

Another possibility, shown in magnetofrictional simulations, is that flux ropes are formed 
and erupt due to slow diffusion of decaying flux from active regions, which interact through 



converging motions on the photosphere ( 


MacKay & van Ballegooijen 


2006 


van Ballegooijen 


& MacKay 


2007 


). Also, 


Chen & Shibata 


(2000 


) simulated the interaction between a flux 



system with a detached flux rope and an emerging flux system that appears to one side of 
the primary fleld. Reconnection between the flux systems leads to a vertical current sheet 
below the detached flux rope, which eventually experiences a fast ejection. 



In this Letter, we consider the 3D extension of the model by Archontis, Hood fc Brady 



(2007) (hereafter called AHB model) in which one (leading) flux tube rises through a stratifled 
atmosphere to create a coronal fleld for the second (following) tube to rise into. The 3D model 
simulates: 1) the intricate reconnection of the two flux systems 2) the multiple eruptions of 
flux ropes, which are initiated by the new emergence of flux and 3) the formation of current 
sheets and hot arcade structures. 



2. Model 

The initial stratiflcation of the stratifled atmosphere and the properties of the magnetic 
flux tubes are similar to the AHB model. The simulations use a 3D version of the Lare shock- 



capturing code (see Arber et al. (2001)). The time-dependent, resistive and compressible 
MHD equations are solved by a Lagrangian remap scheme, including uniform gravity. The 
energy equation is adiabatic and includes viscous and Joule dissipation terms. The resistivity 
is enhanced at sites where the current density exceeds a certain value in a smilar manner 
to the AHB model. At the beginning of the experiment two flux tubes are located 2.9Mm 
and 2.2Mm below the photosphere, with axes parallel to the y direction. The fleldlines are 
uniformly twisted around the main axis of each tube. The initial twist is such that the tubes 
are stable to the kink instability. Both tubes are left-handed and can reconnect when they 
come into contact. The rising motion of the tubes is triggered by a density excess along their 
central axes. The density deflcit is reduced away from = following a Gaussian proflle. 
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similar to the experiment by Archontis et al. (2005). The length scale of the buoyant part 



of the tubes is defined by the parameter A, where A = 20 in our experiments. This deficit 
makes the tubes less dense and more buoyant in the middle part of their axes. As they rise 
through the solar interior, they form f2-shaped loops due to magnetic buoyancy. 

All variables are made dimensionless against their photospheric values, namely: pres- 
sure, P = 1.4 10^ erg cm~^; density, p = 3 10~^ g cm~^; temperature, T = 5.6 10^ K; scale 
height, H = 170 km; time, t = 25 sec; velocity, V = 6.8 km sec~^ and magnetic field, 
B = 1.3 10'^ Gauss. The dimensionless coordinates are (—100 < x < 100), (—100 < y < 100) 
and (— 10 < 2; < 140). The stratified atmosphere includes the top of the solar interior from 
z = —10 to 22, two isothermal layers [z = 22 to 32 for the photosphere and z = 42 to 140 
for the corona) and a transition region with a steep temperature gradient joining them. The 
size of the grid used in the numerical experiments is (256, 256, 320) points in the {x, y, z) 
directions, with x being the transverse coordinate and z the vertical coordinate. 



3. Results and Discussion 

The initial evolution of the system is similar to the evolution in the AHB model. When 
the leading tube reaches the photosphere a bipolar region is formed with a north-south 
orientation due to the initial strong twist of the fieldlines. Eventually, the bipolar region 
moves towards an east- west direction as more internal magnetic layers rise to the photosphere 
and the inclination of the anchored flanks of the tube becomes more vertical. After t = 35, 
the leading tube rises into the non-magnetized corona where it expands and creates an 
ambient field for the following tube to rise into. The direction of the fieldlines and the field 
strength vary away from the center of the expanding volume of the leading system and, 
thus, the newly formed ambient field is more complex than in the previous simulations of 
the interaction between an emerging field and a uniform coronal field. 

As the newly emerged bipolar magnetic field is moving in opposite directions along the 
neutral line, shearing of the field occurs so that the magnetic field lines lose their strong 
azimuthal nature and run nearly parallel to the neutral line. Observations have shown that 
most bipolar sheared magnetic fields in active regions can produce ejective explosions such as 

. Also, in our simulations the twisted and sheared magnetic 
field adopts an overall sigmoidal shape with two oppositely curved elbow-like regions at the 
two (east and west) ends of the neutral line (as illustrated in Fig. [l| left panel). The neutral 
line is stretched along the middle of the sigmoidal structure while the elbows are found on 
either side. During the evolution of the system, the curved region of the elbows consist of 
expanding loops of fieldlines that rise in opposite (north-east and south-west) directions. 



flares and CMEs (Falconer|[2001 
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The parts of the elbows close to the middle of the sigmoid, shear past each other and are 
ready to reconnect if they come into contact. This magnetic field configuration with the 



elbow-like structures has been reported in observations of solar eruptive events by Moore et 



al. (2001). 



When the second flux system emerges at t ~ 40, a thin, curved, sheet-like current 
structure (Fig. [l| middle panel) is formed just above the top layers of the following tube and 
parallel to the interface of the magnetic flux systems, along the y-direction. An interesting 
result is that already at this early stage of the interaction, asymmetry is introduced into the 
system (Fig. [l| left panel) by the 3D nature of the fan-like expansion of the field. On the 
west side of the emerging regions, the fieldlines of the rising systems make a contact angle 
which is larger than on the east side, leading eventually to a higher current concentration 
and more efficient reconnection at the west side of their interface. 

In a similar manner to the 2.5D AHB model, soon after the build up of the interface 
current concentration plasmoid-like structures develop inside the sheet. Due to the asym- 
metry, they are formed firstly towards the west side of the current sheet. The cool and 
dense plasmoids grow and are ejected upwards and sideways, as in the AHB model, allowing 
reconnection in the diffusion regime to occur at a faster rate. 

Two new flux systems are formed by the reconnected fieldlines (Fig. [l| middle panel). 
At the upper end of the current sheet, reconnection forms an overlying field that joins the 
two emerging systems externally, mainly from the north-west side of the leading system 
to the south-east side of the following system while it spreads out mostly above the west 
sides of the emerging systems. At the start of the reconnection, the overlying field links the 
west-side flank of the leading tube with the east-side flank of the following tube. At the 
lower end of the current sheet, reconnection forms a series of magnetic loops in between 
the two emerging systems. The configuration of the low-lying magnetic loops is an arcade- 
like structure joining the leading and following tubes and is reminiscent of post-flare loops. 
The arcade is stretched and curved along the y-direction and grows in size as reconnection 
proceeds and more fieldlines are added to the top of the arcade. The temperature here could 
be as high as 2 MK. The arcade brightens when high velocity reconnection outflows, which 
carry hot material ejected from the current sheet, compress the plasma at the upper part of 
the arcade and increase the temperature. 

During this time, due to the emergence of the leading tube the inner fieldlines are sheared 
horizontally and stretched vertically. Hence, oppositely directed fieldlines, which are rooted 
to the sigmoidal neutral line, are forced into contact. This internal reconnection of fieldlines 



is similar to the 'tether-cutting' mechanism (Moore & Roumeliotis 1992). In experiments 



with a single emerging tube, with the same initial properties, we find that internal 'tether- 
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cutting' reconnection occurs at t ^ 180. In the experiment with two emerging bipoles it 
occurs much earher at t ^ 135. 

In the model with the two tubes, reconnection occurs first (externally) at X-type points 
at the ends of the interface current sheet. At the lower end of the current sheet, the hot 
reconnected fieldlines form an arcade structure. The magnetic pressure inside this new 
arcade increases as reconnection proceeds and eventually exceeds the magnetic pressure of 
the leading flux system. Thus, the arcade expands sideways and pushes the south side of the 
leading system towards its neutral line. It is the combination of shearing along the neutral 
line, stretching due to expansion and the increasing magnetic pressure force of the arcade 
structure that leads to internal reconnection of the fieldlines just above the neutral line of 
the leading emerging system. 

As the two emerging flux tubes continue to rise and interact, internal reconnection 
occurs at many positions along the neutral line of the leading magnetic flux system. Due 
to internal reconnection, flux-rope-like structures are formed inside the magnetized volume 
of the leading tube. Eventually, the ropes rise and erupt into the outer atmosphere. The 
eruption is due to the vigorous sideway pressure by the fast expanding arcade and the release 
of the magnetic tension force of the overlying field, after the ejection of the plasmoids out of 
the interface current sheet. 

A striking result is that the eruptions produced in our simulations do not appear si- 
multaneously. Also, multiple eruptions can be triggered from the same location inside the 
expanding volume of the leading system. As an example, the first formation and eruption of 
a flux rope occurs on the south-west part of the neutral line. This part of the leading tube 
is closer, due to the asymmetry mentioned above, to the following tube. Thus, the dynam- 
ical rearrangment of the magnetic field due to reconnection occurs faster at this region of 
the leading emerging tube. In general, the eruptions occured in our experiment follow the 
direction of the overal 3D expansion of the emerging systems, which for the south-west side 
is preferentially the direction towards the following tube and for the north-east side is the 
opposite direction. Thus, the eruption of the first flux rope is directed towards the south side 
(Fig. [l| right panel). It is a confined eruption, being trapped between the envelope field of 
the leading tube and the expanding field of the upcoming tube that halts the further rise of 
the flux rope. Other eruptions, which occur at later stages of the evolution of the emergence, 
are more dynamic and seem to fully escape from the corona, in a CME-like manner. These 
eruptions are mainly directed towards the north side of the leading magnetic flux system. 

The time evolution of the flux ropes formed inside the leading region, can be followed 
by locating the center (or 0-point) of the ropes at various x z cross-sections (different values 
of y) as functions of time. The height-time relation for 6 values of y is shown in the top 
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panel of Fig. [2j The bottom panel is the projection of their motion onto an x z slice. The 
confined eruption of the first flux rope is shown by the lines at ?/ = — 20 and y = —10. The 
center of the rope for y = —20 is lower than at y = —10 at all times. 

There are two phases during the eruption of the ropes. In the first, the flux ropes 
are moving with moderate speeds, while in the second they are accelerated towards the 
outer atmosphere. These curves are qualitatively similar to the time profiles of filament 



heights shown in the observations by Sterling & Moore (2004, 2005). We also estimate the 



temperature and current density in the area below the central part of the erupting flux ropes 
where a vertical current sheet is formed due to internal reconnection. We find that there is 
a slow increase in the first phase, which is followed by an intensity increase over the period 
of the fast eruption. The peak of the temperature and current density occurs when the fast 
rise is well underway. Since both quantities follow the same time evolution, it is likely that 
the reason for the change in the rise phase of the eruption is that reconnection is slower at 
the beginning but becomes faster, and operates with a higher growth rate, as the ropes rise. 
These results are also consistent with the observations mentioned above in which soft X-ray 
emission is found to occur during the slow-rise phase of a filament's eruption and hard X-ray 
bursts occured after the start of the fast eruption. 

As shown above, the first formation of a flux rope and its eruption occurs in the south- 
west part of the leading tube. Figure [2] (bottom panel) shows that a secondary formation 
and eruption occurs at the same location (y=-10, dashed-triple dotted line), starting at 
t=145. Thus, at least two eruptions are triggered from the same location within the leading 
emerging system until the end of the simulation. Further eruptions may be triggered inside 
the expanding magnetized volumes since the system of the interacting bipoles does not 
reach an equilibrium but evolves dynamically. As the flux ropes erupt, they grow in radial 
extent as additional fieldlines reconnect at the vertical current sheet beneath them. The 
magnetized plasma inside the erupting ropes is a mixture of cold and hot plasma. The cold 
and dense plasma comes from the low atmosphere and is carried upwards by the tether- 
cutting reconnection forming concave upwards fieldlines. The hot plasma (up to 2 MK) is 
emitted inside the envelope of the rope by the reconnection jet, which blasts from the upper 
end of the vertical current sheet below the rope. The velocity at the center of the rope is 
^10-15Km/sec. However, the maximum upward velocity inside the complete volume of the 
rope reaches values up to 150 Km/sec, mainly due to the emission of the reconnection jet. 
The rising motion of the flux rope can result in a CME-like eruption. 



Figure |3](left panel) shows the temperature distribution in a vertical plane and magnetic 
fieldlines, which are traced from near the axis of the flux rope. The cold plasma is mostly 
concentrated near the axis of the flux rope while the hot plasma is overlying the axis on the 
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right side of the eruption. The fieldhnes near the axis of the rope come from both of the 
emerging flux systems. More precisely, fieldhnes from the following tube joins the leading 
tube, through the arcade structure between them before forming the erupting flux rope. 
These fieldhnes undergo multiple reconnection events along their lengths and, thus, may 
end up in either of the flanks at the east-side boundary. The middle panel in Fig. [3] shows 
isosurfaces of current density, which are overplotted on the same set of fieldhnes. A striking 
result is that long and twisted sheets are found close to the core of the erupting flux rope. In 
addition, arch-like current sheets that are less twisted and have different orientation to each 
other, are formed in the region between the emerging systems (see 'A' and 'B' in the figure). 
These structures, with enhanced current density, may be quasi seperatrix layers. The right 
panel in Figure [3] is a blow-up of the region at the vertical current sheet below the erupting 
flux rope shown in the left panel. The fieldhnes shown in the figure have been traced from the 
lower end of the vertical sheet. These are reconnected fieldhnes that surround the leading 
flux tube's original axis. Thus, the field below the vertical current sheet forms a tightly 
wound flux tube with its axis lying parallel to the neutral line of the leading bipolar region. 
The temperature rises up to 0.2 MK at the bottom of the current sheet, as hot reconnected 
plasma moves downward from the sheet. These fieldhnes are also linked: a) to fieldhnes that 
envelope the erupting flux rope from the north side (two fieldhnes on left side of figure) and 
b) to fieldhnes that form the hot arcade structure between the emerging flux systems. 

Our 3D simulations show that the interaction of emerging flux tubes may lead to multiple 
eruptions of flux ropes. A network of twisted current sheets, with an overall S-shaped 
configuraton, and hot cusp-like, arcade structures are also formed during the evolution of 
the system. 

Computational time on the Linux clusters in St. Andrews (STFC and SRIF funded) 
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Fig. 1. — Left: Visualization of the full magnetic field vector (arrows) at photospheric heights, at 
t = 130. Also shown is fieldlines of the emerging tubes and isosurfaces of current density. The 
region around the asymmetry is underlined by the inset in the figure. Middle :Fieldlme topology 
around the current sheet at the interface at t = 135. The darkest the color the weakest the magnetic 
field is along the fieldlines. Right: Confined eruption of a flux rope at t = 140. Superimposed is 
the magnetic field vector (arrows) in a vertical plane at y = —22. Sites of internal reconnection are 
shown in the inset rectangle of the panel. 
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Fig. 2. — Top: Height-time relation of the centers of the flux ropes. Bottom: Projection of the 
rising motion of the centers of the flux ropes on a horizontal plane. 




Fig. 3. — Left: Ejective eruption of a flux rope at t = 175. The temperature contours are on 
the plane y = 20. Middle: The same as in the left panel with isosurfaces of current density 
superimposed. Right: 3D visualization of fieldlines in the region below the erupting flux rope. 
Superimposed on the temperature contours is the velocity field (arrows) on y = 20. 



